Role of electronic correlations in Ga 
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Abstract 

An extended around mean field (AMF) functional for less localized p electrons is developed to 
quantify the influence of electronic correlations in a-Ga. Both the local density approximation 
(LDA) and generalized gradient approximation (GGA) are known to mispredict the Ga positional 
parameters. The extended AMF functional together with an onsite Coulomb interaction of U e s = 
1.1 eV, as obtained from constraint LDA calculations, reduces the deviations by about 20%. The 
symmetry lowering coming along with the electronic correlations turns out to be in line with the 
Ga phase diagram. 
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Electronic correlations beyond the local density and generalized gradient approximation 
(LDA/GGA) to density functional theory (DFT) are usually considered to be exclusive to 
partially filled d and / states [l^j]. The localized nature of these orbitals results in a non- 
negligible onsite Coulomb interaction U 5|-|9|]. The importance of correlations in partially 
filled p orbitals recently has been pointed out in the context of u d° magnetism" , where the 
magnetic order arises from the p states instead of the conventional d and / states lOj, 111 ]- 



First-principles calculations have shown that correlations in open p shells exist not only for 
first-row elements, like N and O (2p), but also for heavier atoms, like Te {hp) [l2,[l^]. Up 
to now, the interest in p orbital electronic correlations was limited to pure and doped ionic 
semiconductors. 

Great interest in a-Ga was triggered by the common belief that it is the only elemental 
solid in which metallicity and covalency coexist jl^| . Showing orthorhombic crystal symme- 



try with space group Cmca, a-Ga is the most stable Ga phase at ambient conditions |15l. |16|. 
Besides the lattice constants a = 4.5102 A, b = 4.5161 A, and c = 7.6448 A, the crystal 
structure is determined by the two internal parameters u = 0.0785 and v = 0.1525, which 
specify the Ga atomic positions in the unit cell 17|. Every Ga atom (Ga re f in Fig. [TJ has 
one nearest neighbor (Gai) and six neighbors in the next three coordination shells (Ga2 to 
Ga4). Partial covalency of a-Ga, indicated by the short Ga re f-Gai bond, highly anisotropic 
electronic and thermal conductivities, and a steep pseudogap around the Fermi level have 




FIG. 1. (Color online) Local environment of a reference atom Ga re f. Gaj (i = 1,2,3,4) denotes 
the atoms in the i-th coordination shell of Ga re f. Notice that there are two Ga2 atoms due to the 
periodicity along the 6-axis. 
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been observed in experiments [l8M22j| and LDA/GGA calculations . The different 



bond lengths between Ga re f and Ga« (i = 1,2,3,4) constitute a special local environment, 
which is crucial to the coexistence of metallicity and covalency in a-Ga. As it is reported 



in Ref. jl7| and confirmed by our calculations, however, the LDA/GGA fails miserably to 
reproduce the experimental Ga positional parameters. It is our aim to evaluate to which 
extent electronic correlations account for this deficiency. 

The standard way to treat electronic correlations in a DFT calculation is to consider an 
onsite parameter U and to embed a Hubbard-like Hamiltonian into the LDA/GGA Kohn- 
Sham equations (LDA/GGA+£7 method) [2]. However, usually the onsite U is applied only 
to the electrons inside the non-overlapping muffin-tin spheres. By the neglection of effects 
from the interstitial electrons, any result will depend on the muffin-tin radius R m t. Major 
inaccuracies are expected for more delocalized states, like the 4p electrons of a-Ga. As a 
consequence, an extended LDA/GGA+f/ functional is required which includes the effects 
of the interstitial charge. Due to the reduced localization of the p electrons, a less strongly 
correlated state is expected. Hence, the around mean field (AMF) approximation is chosen 



for modeling a-Ga 



24]. 



Based on the fact that the LDA corresponds to the homogeneous solution of the mean- 
field Hartree-Fock equations with equal occupancy of all sub-orbitals with the same spin, 
the AMF energy functional is obtained by supplementing the LDA/GGA functional by the 



additional term 



£ AMF = -^£Tr( n(7 -rU) 2 . (1) 



In this relation U e g = U — J is the effective interaction, where U and J denote the onsite 
Coulomb and exchange interactions, respectively, and n a the density matrix for spin a. In 
addition, n a = Tr(n CT )/(2£ + 1) is the average occupation number of the sub-orbitals with 
spin a (and orbital quantum number £) and I the identity matrix. The double counting is 
already corrected in Eq. (jTJ). 

We assume a linear dependence = X a n™ f ' with an orbital dependent parameter \ a . 
Because n CT and n™* can be obtained from the Mulliken method and the DFT calculations, 
respectively, \ a = n CT (n° lt ) _1 can be evaluated. For the Mulliken analysis we will apply the 



Gaussian code 



25) . In general, A CT will depend on the choice of the basis set. However, we 



have checked a series of basis sets to ensure that the dependence is sufficiently weak. To fix 



the dependence on R mt , we choose Tr\ a j (2/ + 1) = 2.67 with R mt = 2.2 a.u. in the following. 
The (orbital dependent) one-electron potential now is given by the functional derivative of 
E AMF with respect to n™*, 



Vr' mt = -t/A(n ff -n ff I). (2) 
We have implemented the extended AMF approach, given by Eqs. ([T]) and ([2]), for the 



WIEN2k code 24J. Like the original, the modified one-electron potential is applied to this 
part of a sub-orbital that is inside the muffin-tin sphere. Yet, there are two improvements: 
First, the potential exerted by a particular sub-orbital arises from the Coulomb repulsion of 
electrons both inside and outside the muffin-tin sphere, n m CT , instead of only n^ a . Second, 
the parameter \ m ,a mimicks the energetical effects which would be obtained if both parts 
of the sub-orbital (inside and outside the muffin-tin sphere) were subject to the potential. 
I.e., the extended AMF functional includes effects of charge not only inside the muffin-tin 



TABLE I. Fully optimized structural parameters (a, b/a, c/a, u and v) of a-Ga obtained from 
standard LDA/GGA calculations. Different XC functionals, including (1) LDA-PW92, (2) GGA- 
SOGGA, (3) GGA-AM05, (4) GGA-PBEsol, (5) GGA-WC, (6) GGA-PBEalpha, (7) GGA-PBE, 
(8) GGA-PW91, (9) GGA-BPW91, and (10) GGA-RPBE, are used Q- The experimental values 
are taken from Ref. [171 ]. 





a (A) 


b/a 


c/a 


u 


V 


SVol (%) 


Exp 


4.5102 


1.695 


1.0013 


0.0785 


0.1525 





(1) 


4.4262 


1.695 


1.0014 


0.0843 


0.1559 


-5.48 


(2) 


4.4378 


1.705 


1.0074 


0.0821 


0.1558 


-3.59 


(3) 


4.4646 


1.703 


1.0064 


0.0815 


0.1561 


-2.03 


(4) 


4.4658 


1.703 


1.0064 


0.0829 


0.1559 


-1.95 


(5) 


4.4816 


1.698 


1.0033 


0.0834 


0.1559 


-1.51 


(6) 


4.5300 


1.698 


1.0033 


0.0825 


0.1561 


1.72 


(7) 


4.5687 


1.698 


1.0033 


0.0803 


0.1563 


4.34 


(8) 


4.5743 


1.695 


1.0013 


0.0805 


0.1564 


4.32 


(9) 


4.5839 


1.695 


1.0013 


0.0826 


0.1562 


4.97 


(10) 


4.6344 


1.695 


1.0013 


0.0808 


0.1566 


8.48 
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spheres but also in the interstitial region. Thus, it eliminates the dependence of the results 
on R mt , making it suitable to study electronic correlations in less localized p orbitals. 

Using standard LDA/GGA, we have fully optimized the structure of a-Ga against all 5 
structural parameters, applying Pulay corrections to the forces. We use in all calculations 
i? mt = 2.2 a.u., -Rmt-Kmax = 8.5, and l max = 10 as well as the same 23 x 23 x 14 /c-mesh. 
The validity of the LDA/GGA is tested for ten exchange- correlation (XC) functionals, see 
Table [H In order to provide a quantitative comparison between the experimental findings 
and our results for the unit cell volume (Vol) we study SVol = (Vol — Vo1exp)/Vo1e xp - The 
a lattice parameter depends strongly on the choice of the XC functional, while both ratios 
b/a and c/a are almost constant. As a consequence, SVol varies within a wide range from 
-5.48% to 8.48%. The LDA-PW92 functional yields the smallest value of a and thus the 
smallest volume, which reflects the common insight that the LDA, in contrast to the GGA, 
underestimates the volume . The overestimation of u and v in all our data indicates 

a systematic failure of the LDA/GGA to describe a-Ga. This may be due to the influence 
of electronic correlations, which we probe in the following. 

We apply our extended AMF functional to the outer 4p orbitals to optimize the crystal 
structure, choosing the GGA-WC XC functional for which SVol is minimal. The computa- 
tional details are the same as in the LDA/GGA calculations. To establish the strength of 
the onsite Coulomb interaction, we employ the constraint LDA method [27] , which leads to 
U e g = 1.1 eV. With this value and the experimental lattice constants, u and v are optimized 
in the extended AMF scheme. The anisotropics of the electronic conductivity (a c :ab.a a ) and 
thermal conductivity (K c :Kb:K a ) are obtained [281], see Table HT1 The extended AMF calcu- 
lations lead to a qualitative improvement of both the structural and transport properties. 
Most remarkably, u is improved by about 20% with respect to the experimental value. The 



TABLE II. Experimental 13, LL8| an d calculated positional parameters and transport anisotropies. 



Exp U eS = eV U eS = 1.1 eV 

u 0.0785 0.0835 0.0824 

v 0.1525 0.1560 0.1556 

a c :a b :cr a 1:3.1:6.7 1:1.8:4.9 1:2.0:5.3 

K c :K b :K a 1:3.0:6.6 1:2.0:4.4 1:1.8:4.7 
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FIG. 2. (Color online) (a)-(e) Total and partial DOS obtained for the GGA-WC (black solid line) 
and extended AMF (red dotted line) methods, (f) Charge density difference An in the (100) plane 
at x = 0.5. The positions of the Ga atoms are indicated by black dots. 



extended AMF calculations hence predicts a shorter Ga re rGa! bond and a slightly higher 
anisotropy of the electronic and thermal conductivities. I.e., the partial covalency of cv-Ga, 
which is underestimated by the GGA-WC, is improved. Still, the total and partial s, p x , p y , 
and p z densities of states (DOS) in Figs. Hfa)-(e) reveal only little changes. 

On the other hand, an enhanced covalency along the c-axis due to the onsite interaction 
is reflected by the charge density. In Fig. [2](f ) we show the charge density difference within 
the (100) plane, at x — 0.5, between the GGA-WC and the extended AMF results: An = 
riAMF — ^gga-wc- An accumulation of charge is found in the region between Ga re f and Gai, 
indicating a stronger Ga re f-Gai bond and, thus, an enhanced covalency. The accumulation 
comes along with charge transfer from the p x and p y sub-orbitals to the p z sub-orbital. The 
p x and p y occupations decrease from 0.176 and 0.187 to 0.174 and 0.185, respectively, while 
the p z occupation grows from 0.204 to 0.208. The principal component of the electric field 
gradient, the direction of which is indicated by the arrow in Fig. [2](f), grows from 4.2 • 10 21 
V/m 2 to 5.3 • 10 21 V/m 2 , whereas the asymmetry parameter decreases from 0.207 to 0.153. 
The experimental value is 5.5 ■ 10 21 V/m 2 with an asymmetry parameter of 0.179 29J. The 
finite U e g of the Ap electrons is related to the anomalous spatial contraction of the valence 
orbitals in a-Ga, which is a consequence of incomplete screening of the nuclei by a relatively 
shallow Ga 3d state 30j. For the same reason, also a localization of the 4s electrons would 
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FIG. 3. (Color online) Relation between asymmetry, characterized by Pasymmj and stability, char- 
acterized by the pressure at which a phase occurs, for various Ga crystalline phases. The data are 



taken from Refs. 



ia 



31 



321 ] . The red dotted line is a guide to the eye. 



be expected. Unlike the 4p electrons, however, the Ueg has no effect on the 4s electrons in 
the AMF approximation 2j. 

By Eq. (CQ), a symmetry lowering due to electronic correlations is expected. Because an 
asymmetric crystal structure gives asymmetric orbital populations, p a symm = {J2 m A( n m,a — 
n a ) / n a } 2 }^ , the electronic correlation energy is related to the asymmetry of the structure by 
the relation E AMF = — U e gfi 2 p 2 symm /2, where n = = fi± in non-magnetic a-Ga. Electronic 
correlations thus stabilize a crystal structure with lower symmetry, due to a smaller _£; AMF . 
This picture is confirmed by the successful prediction of the quadrupolar lattice distortion 
in the perovskite compound KCuF 3 [3], where the introduction of the U parameter creates 
a "mexican-hat" shaped energy surface. This shifts the energy minimum towards a lattice 
with lower symmetry. The symmetry lowering due to electronic correlations is also reflected 
by the orbital polarization of the extended AMF potential functional in Eq. ([2]). If n mjCr is 
larger (smaller) than n a , it further increases (decreases) due to a smaller (larger) V^f F ' mt . 
The orbital polarization manifests in the aforementioned occupation numbers of a-Ga. 

If electronic correlations are present in a-Ga, the symmetry effect should be seen in the 
Ga phase diagram. The Ga-IV phase, which has a face-centered cubic structure, occurs only 
at a pressure above 120 GPa [311 ] . In contrast, the a-Ga phase, which has an orthorhombic 



structure with the lowest symmetry, is the most stable phase at ambient conditions |15l. |16|. 
For quantitative description, we use the pressure at which a phase occurs as measure of its 
stability, i.e., higher/lower pressure corresponds to lower/higher stability. In Fig. |3]we show 
Pasymm as a function of pressure. The results confirm that structures with lower symmetry 
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give rise to more stable phases. 

In conclusion, we have extended the AMF functional to describe correlation effects of less 
localized p electrons. The partial covalency of a-Ga, which is underestimated in LDA/GGA 
calculations, is improved by our approach. This shows that finite electronic correlations are 
present in the partially filled 4p orbitals. They are reflected by the crystal symmetry of the 
different phases in the Ga phase diagram. 
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